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Technical  Section 


Technical  Objectives 

The  objectives  of  the  grant  were  to  provide  a  systematic  study  to  fill  the  gap  between  existing  research  on 
low  Reynolds  number  turbulent  flows  to  the  kinds  of  turbulent  flows  encountered  on  full-scale  vehicles. 
We  report  specifically  on  (1)  the  behavior  of  wakes  at  high  Reynolds  numbers  with  increasing 
complexity,  varying  from  axisymmetric  wakes,  to  wakes  typical  of  appended  undersea  platforms;  and  (2) 
the  effects  of  roughness  at  high  Reynolds  numbers.  We  have  gained  a  better  understanding  of  complex 
flow  interactions  in  wakes  typical  of  naval  platform  flow  fields,  especially  wakes  with  concentrated 
regions  of  streamwise  vorticity.  We  have  also  made  progress  in  the  understanding  of  high  Reynolds 
number  flows  over  rough  surfaces,  especially  pipes  and  flat  plates  with  roughness  that  relate  to  marine 
surfaces.  We  expect  these  studies  to  lead  to  improved  flow  prediction  and  improved  flow  control.  The 
work  was  performed  in  two  unique  facilities:  the  Superpipe  and  the  High  Reynolds  number  Test  Facility 
(HRTF)  that  can  obtain  very  high  Reynolds  numbers  on  a  laboratory  scale  using  compressed  air  as  the 
working  fluid. 


Technical  Approach 

Our  studies  are  performed  in  two  facilities  constructed  at  Princeton  under  ONR  funding.  The  Supeipipe 
facility  enables  very  accurate  measurements  in  fully  developed  turbulent  pipe  flow  across  a  wide  range  of 
Reynolds  numbers,  from  31  x  103  to  35  x  106.  High  Reynolds  numbers  are  achieved  at  a  moderate  cost 
by  using  compressed  air  at  ambient  temperatures  as  the  working  fluid,  thereby  decreasing  the  kinematic 
viscosity  by  over  two  orders  of  magnitude  as  compared  to  air  at  STP.  The  maximum  static  pressure  is 
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200  atm.  The  test  pipe  has  a  nominal  diameter  of  129  mm,  and  a  length  of  202  D.  A  diagram  of  the 
facility  is  shown  in  Figure  1. 
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Figure  1:  Princeton/D  ARP  A/ONR  Superpipe  apparatus. 

A  wind  tunnel  to  achieve  high  Reynolds  numbers  called  the  Princeton/ONR  High  Reynolds  Number 
Testing  Facility  (HRTF)  is  also  used  to  attain  high  Reynolds  number.  Like  the  Superpipe,  it  uses  air  at 
pressures  up  to  3,000psi  as  the  working  fluid  (Figure  2).  The  primary  purpose  of  the  facility  is  to  study 
the  hydrodynamic  forces,  moments  and  flow-fields  produced  by  advanced  naval  platforms  up  to  length 
Reynolds  numbers  of  150  x  106.  There  are  two  working  sections:  each  is  8ft  long  with  an  internal 
diameter  of  18in.  The  facility  will  be  equipped  with  a  Magnetic  Suspension  Balance  System  (MSBS)  to 
allow  measurements  free  of  the  interference  produced  by  the  support  systems  usually  employed  in  these 
applications.  The  MSBS  is  being  designed,  constructed  and  tested  under  the  supervision  of  Colin  Britcher 
ofODU. 
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Progress 


1.  Honed  rough  pipe  flow  mean  flow  experiments 

Mean  flow  measurements  are  reported  by  Shockling  et  al.  (2006)  and  Allen  et  al.  (2006).  These 
experiments  have  shown  that  honed  roughness  has  an  inflectional  transitional  roughness  behavior,  similar 
to  sand-grain  roughness,  and  unlike  the  monotonic  behavior  expected  according  to  Colebrook’s 
transitional  roughness  function,  on  which  the  Moody  Diagram  is  based  (see  Figure  3).  In  combination 
with  the  work  of  Schultz  and  Flack  at  the  Naval  Academy,  we  tentatively  propose  that  inflectional 
roughness  is  much  more  common  than  previously  known,  and  may  well  form  the  correct  model  for  ship 
roughness.  Measurements  of  the  honed  surface  roughness  function  show  that  roughness  begins  to  have 
an  effect  on  the  velocity  profile  when  the  equivalent  sandgrain  roughness  exceeds  about  3.5  (Figure  4). 
Outer  flow  similarity  as  proposed  by  Townsend  was  confirmed  at  all  Reynolds  numbers  (Figure  5). 


Figure  3:  Friction  factor  for  the  present  surface,  compared  with  the  rough-wall  relations  of  Colebrook 
(1939)  for  the  same  equivalent  sandgrain  roughness,  the  smooth- wall  relation  of  McKeon  et  al.  (2005), 
and  the  results  for  the  smallest  sandgrain  roughness  used  by  Nikuradse  (1933). 


Figure  4:  Left:  Velocity  profiles  over  the  full  Reynolds  number  range.  Right: . ,  Hama  roughness 

function. - ,  Colebrook  roughness  function  with  ks  =  7.4  x  10'6  m. 
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2.  Honed  rough  pipe  flow  turbulence  experiments 

The  hot-wire  measurements  of  the  streamwise  component  of  turbulence  have  been  completed  and  are  in 
press  (Kunkel  et  al.  2007).  Streamwise  turbulence  intensities,  higher  order  moments,  and  spectra  were 
compared  with  the  corresponding  results  from  a  previous  smooth  wall  pipe  flow  study  in  the  same  facility 
(Figures  6,  7  and  8).  The  hydrodynamically  transitionally  rough-wall  turbulence  intensities  and  higher- 
order  moments  agreed  well  with  the  smooth-wall  results  indicating  that  the  roughness  does  not  affect  the 
time  averaged  statistics  in  the  outer-region  of  the  flow.  The  rough-wall  one-dimensional  streamwise 
spectra  in  the  logarithmic  region  are  also  found  to  agree  with  the  smooth-wall  results  indicating  that  the 
roughness  does  not  affect  the  distribution  of  energy  in  the  flow  and  implies  that  the  it  does  not  change  the 
structure  of  the  flow. 


Figure  6:  Streamwise  turbulence  intensities.  Solid  symbols  are  rough  pipe  data  and  hollow  symbols  are 
smooth  pipe  data  reproduced  from  Morrison  et  al.  (2004).  (a)  Outer  layer  scaling;  (b)  Inner  layer  scaling. 

Overall  these  data  support  Townsend's  Reynolds  number  similarity  hypothesis  in  that  the  energy- 
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containing  motions  are  independent  of  roughness  aside  from  the  increased  wall  shear  stress.  The  current 
measurements,  while  unique  and  insightful,  are  complicated  by  probe  resolution  problems  at  low  wall- 
normal  positions  and  high  wave  numbers.  Future  work  will  attempted  to  resolve  these  issues  (see  Kunkel 
et  al.  2007)  as  well  as  measure  the  wall-normal  velocity  component  of  the  flow.  While  other  authors  have 
found  slight  differences  in  the  streamwise  turbulence  statistics,  the  wall-normal  statistics  were  shown  to 
be  the  most  sensitive  to  wall  roughness. 


Figure  7:  Smooth-  and  rough- wall  pipe  (a)  skewness,  (b)  flatness.  Symbols  as  in  Figure  6. 


Figure  8:  Rough  and  smooth  pipe  streamwise  velocity  spectra  at  ReD,  R"  =  5.5  x  106,  1.0  x  105.  Solid 
symbols  are  rough  pipe  data  and  hollow  symbols  are  smooth  pipe  data  reproduced  from  Morrison  et  al. 
(2004).  (a)  Outer  layer  scaling;  (b)  Inner  layer  scaling. 


3.  Commercial  steel  rough  pipe  flow  experiments 

The  honed  pipe  described  above  was  replaced  with  a  commercial,  extruded  steel  pipe.  The  surface 
roughness  was  as  delivered  by  the  supplier,  with  k^  =  5  pm.  Results  on  friction  factor  and  velocity 
profiles  are  shown  in  Figures  9  to  11.  These  are  the  first  results  ever  obtained  in  commercial  steel  pipe 
under  laboratory  conditions.  They  show  that  commercial  steel  pipe  has  a  friction  factor  behavior  that  lies 
between  the  Colebrook  type  roughness,  and  sand-grain  type  roughness,  in  that  the  departure  from  the 
smooth  pipe  is  rather  sudden,  with  little  or  no  evidence  for  the  “inflectional”  behavior  characteristic  of 
sand-grain  roughness  (and  honed  surface  roughness).  In  addition,  the  equivalent  sand-grain  roughness  of 
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commercial  steel  pipe  is  closer  to  2.0kmis,  rather  than  the  commonly  accepted  value  of  3.0krms.  The 
velocity  profiles  in  the  rough  regime  show  the  expected  departure  below  the  logarithmic  law,  and  they 
also  follow  Townsend’s  outer-layer  similarity  (as  did  the  honed  pipe).  The  work  has  been  submitted  for 
publication  (Langelandsvik  et  al.  2007)). 


Superpipe,  natural  rough  steel  pipe 


R«H 

Figure  9:  Friction  factor  for  the  commercial  steel  pipe  surface,  compared  with  the  rough-wall  relations  of 
Colebrook  (1939)  for  the  same  equivalent  sand-grain  roughness,  and  the  smooth-wall  relation  of  McKeon 
etal.  (2005). 


Velocity  profile,  inner  scaling,  low  Re 
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Figure  10:  Velocity  profiles  in  commercial  steel  pipe  for  Reynolds  numbers  up  to  220  x  10'3. 


6 


Figure  1 1 :  Velocity  profiles  in  commercial  steel  pipe  for  Reynolds  numbers  up  to  3  x  10'6. 


4.  Boundary  layer  experiments 

Additional  tests  have  been  performed  on  a  flat  plate  boundary  layer  to  establish  the  use  of  tunnel  inserts  to 
produce  a  zero  pressure  gradient  environment  (Figures  12  and  13).  These  experiments  are  preparatory  to 
the  HRTF  experiments  on  the  same  plate.  The  mean  velocity  (Figure  13),  turbulence  intensity  (Figure 
14),  and  spectra  (Figure  15)  all  demonstrate  that  the  experimental  arrangements  used  in  the  HRTF  to 
study  flat  plate  boundary  layers  are  valid  to  an  extraordinary  degree. 


Figure  12:  Boundary  layer  experiment  on  the  flat  plate.  “Without  lid”  means  in  the  empty  tunnel  (left). 
“With  lid”  means  results  on  the  flat  plate  with  a  semi-circular  cover  that  includes  an  insert  to  produce  a 
zero  pressure  gradient  (right). 
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Figure  13:  Boundary  layer  results  on  the  flat  plate.  “Without  lid”  means  in  the  empty  tunnel.  “With  lid” 
means  results  on  the  flat  plate  with  a  semi-circular  cover  that  includes  an  insert  to  produce  a  zero  pressure 
gradient. 


Turbulence  Intensity,  Inner  Scaling 
Re_theta  ~  6,000 


Figure  14:  Turbulence  intensity  in  boundary  layer  on  the  flat  plate.  “Rectangular  cross-section”  means  in 
the  empty  tunnel.  “Semi-circular  cross-section”  means  results  on  the  flat  plate  with  a  semi-circular  cover 
that  includes  an  insert  to  produce  a  zero  pressure  gradient. 
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Premultiplied,  Inner  log  region,  inner  scaling 
Re~5,000 


k1*y 


Figure  15:  Spectra  in  turbulent  boundary  layer  on  the  flat  plate.  “Rectangular  cross-section”  means  in  the 
empty  tunnel.  “Semi-circular  cross-section”  means  results  on  the  flat  plate  with  a  semi-circular  cover  that 
includes  an  insert  to  produce  a  zero  pressure  gradient. 


5.  SUBOFF  wake  experiments  for  the  HRTF 

The  model  has  been  completely  pressure-tapped  for  static  pressure  measurements,  and  equipped  with  a 
Scanivalve  system.  The  wake  traversing  system  (2  degrees  of  translation,  2  degrees  of  rotation)  has  been 
completed  and  installed  in  the  HRTF  working  section  (Figures  16,  17  and  18).  Additional  screen  and 
mesh  have  been  installed  to  improve  flow  quality  and  reduce  freestream  turbulence  level),  as  well  as  the 
model  yaw  motion  system. 


Figure  16:  Traversing  system  as  installed  in  the  HRTF  downstream  of  SUBOFF  model. 
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Figure  18:  HRTF  control,  power,  and  pressure  transducer  port.  Left:  disassembled.  Right:  assembled. 


Measurements  were  made  in  the  wake  of  a  DARPA  SUBOFF  model  over  a  large  range  of  Reynolds 
numbers  based  on  model  length,  from  1.1  x  106  to  25  x  1 06,  on  the  centerline  of  the  wake  for  locations  3, 
6,  9,  12,  and  15  diameters  downstream  from  the  tail.  The  model  was  axisymmetric  without  appendages 
(fins)  and  it  was  supported  by  a  streamlined  support.  That  is  a  semi-infinite  extension  of  the  sail.  The 
wake  experimental  results,  obtained  using  Particle  Image  Velocimetry  (PIV)  and  crossed  hot-wires, 
provided  qualitative  and  quantitative  insight  into  the  flow  field  created  by  an  undersea  platform.  In 
addition,  the  pressure  was  measured  at  45  different  locations  along  the  model  for  three  different  Reynolds 
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numbers,  1.1  x  106,  12  x  106  and  25  x  106.  Also,  PIV  measurements  were  conducted  in  the  wake  of  the 
sail  attached  to  a  DARPA  SUBOFF  model  at  9.36  x  105.  Four  different  yaw  angles,  between  6  and  17 
degrees  were  investigated  yielding  insights  into  the  behavior  of  the  junction/hull  and  sail  tip  vortices. 

Mean  flow  and  turbulence  data  obtained  using  crossed- wire  anemometry  for  Re  =  25  x  106  are  shown  in 
figures  19,  20  and  21.  For  all  Reynolds  numbers  studied,  the  mean  velocity  distribution  attains  self¬ 
similarity  at  distances  between  3  and  6  diameters  downstream  for  the  side  where  the  support  is  not 
located,  and  follows  an  exponential  function  as  expected  from  similarity  arguments.  In  contrast,  the  mean 
velocity  distribution  for  the  support  side  does  not  attain  self-similarity,  and  displays  significant  effects  of 
the  support  wake  and  support/body  junction  flows.  In  addition,  none  of  the  Reynolds  stress  distributions 
of  the  flow  attain  self  similarity.  For  the  higher  Reynolds  numbers  studied  the  presence  of  the  support 
introduces  an  asymmetry  into  the  wake  that  results  in  the  overall  decrease  of  radial  and  axial  turbulence 
intensities  for  the  support  side.  Also  for  the  higher  Reynolds  numbers,  the  axial  turbulence  intensity 
distribution  is  relatively  symmetric  while  maintaining  a  maximum  on  the  side  without  the  support.  The 
opposite  is  observed  for  Re  =  1 .1  x  106.  Furthermore,  the  values  of  the  turbulent  shear  intensity,  u'v',  at 
Re  =  1.1  x  106  are  higher  on  the  support  side,  r/D  <  0,  for  each  x/D  location  measured.  This  trend  is 
reversed  for  the  two  higher  Reynolds  number  cases  with  the  only  exception  of  x/D  =  3 .  Also,  the 
coefficient  of  pressure  distribution  along  the  top  meridian  line  of  the  model,  r/D  >  0,  is  generally  lower 
for  Re  =  1.1  x  106  than  that  for  Re  =  12  x  106  and  25  x  106,  which  seem  to  have  collapsed.  The  sail  wake 
experiments  demonstrated  the  significance  of  the  sail  tip  vortex  when  the  model  is  at  a  nonzero  yaw 
angle.  As  the  yaw  angle  is  increased  from  6  to  12  degrees  the  circulation  of  the  sail  tip  vortex  increases. 
As  the  yaw  angle  is  further  increased  the  boundary  layer  separates  with  an  overall  drop  in  circulation.  A 
similar  phenomenon  is  observed  for  the  junction  vortex  with  the  exception  that  when  the  yaw  angle  is 
further  increased  to  17  degrees  the  circulation  continues  to  increase  at  a  slower  rate.  Also,  the  circulation 
values  for  the  sail  tip  vortex  are  about  an  order  of  magnitude  larger  than  those  of  the  junction  vortex. 


Figure  19:  Mean  flow  results  in  similarity  coordinates  for  Re  =  25  x  106. 

The  effects  of  the  support  on  the  wake  development  are  similar  to  the  effects  introduced  by  the  sail  on  a 
model  wake  (except  for  the  absence  of  the  tip  flow).  The  presence  of  the  support  affects  the  flow 
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differently  for  different  Reynolds  numbers  emphasizing  the  importance  of  high  Reynolds  number  studies 
to  better  understand  naval  platform  flows. 


-4  -2  0  2  4 


r/ la 

Figure  20:  Streamwise  turbulence  results  in  similarity  coordinates  for  Re  =  25  x  106. 
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6.  Magnetic  Suspension  Balance  System  (MSBS) 

Work  by  Colin  Britcher  at  ODU:  2  degree-of-freedom  levitation  has  been  achieved  inside  the  test  section; 
5  degree-of-freedom  levitation  inside  is  imminent  (Figure  22).  Colin  Britcher  will  report  this  progress 
separately. 


Figure  22:  MSBS  system  as  currently  configures  (presently  located  at  NASA  LaRC). 


Technology  Transfer 

The  Norwegian  company  Gassco  has  supported  a  graduate  student  (Leif  Langelandsvik)  for  6  months  on 
studies  of  roughness  in  the  Princeton  Superpipe  apparatus,  in  partial  support  of  the  ONR  program.  Gassco 
is  continuing  discussions  on  future  collaborations. 
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